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Abstract 

In support of application of the DARPA-AIM methodology to the 
accelerated hybrid thermal process optimization of 3 rd generation 
aeroturbine disc alloys with quantified uncertainty, equilibrium 
and diffusion couple experiments have identified available 
fundamental thermodynamic and mobility databases of sufficient 
accuracy. Using coherent interfacial energies quantified by 
Single-Sensor DTA nucleation undercooling measurements, 
PrecipiCalc ™ simulations of nonisothermal precipitation in both 
supersolvus and subsolvus treated samples show good agreement 
with measured y’ particle sizes and compositions. Observed long- 
term isothermal coarsening behavior defines requirements for 
further refinement of elastic misfit energy and treatment of the 
parallel evolution of incoherent precipitation at grain boundaries. 

Introduction 

An important ongoing application of predictive science-based 
computational materials engineering is the accelerated hybrid 
process optimization of dual micro structure aeroturbine discs 
fabricated from 3 rd generation high-refractory Ni-base 
superalloys. In support of higher fidelity application of the 
DARPA-AIM methodology previously demonstrated on the 
IN 100 and R88DT disc alloys [1,2], a collaborative model and 
database validation study is supported by the NASA Aviation 
Safety Program to better quantify uncertainty and improve 
prediction accuracy of the QuesTek PrecipiCalc 
multiphase/multicomponent precipitation simulation code 
developed as the primary engine of the AIM methodology. 


Table 1 Nominal compositions (balance Ni) of the four 3 rd 
generation disc alloys studied under this program. 



Four 3 rd -generation nickel-based disk superalloys, including ME3 
(also called Renel04), LSHR (Low-Solvus, High-Refractory 
alloys developed by NASA), Alloy 10 (developed by Honeywell), 
and RR1000 (developed by Rolls-Royce) are being studied under 
this program, and their nominal compositions are listed in Table 1. 
In the following sections, key results from the 1 st year of this 
NASA/QuesTek collaboration are provided. 


PrecipiCalc Calibration and Validation Protocol 

An initial objective of this NASA program is to develop a 
standard PrecipiCalc calibration and validation protocol for 
intragranular y’ precipitation. The protocol needs to employ 
independent experimental measurements to decouple and/or 
minimize the cross interaction between model parameters, 
allowing the determination of the model parameters with high 
fidelity and minimum overfitting. The protocol must sequentially 
address the foundational databases and model parameters of the 
PrecipiCalc method, which include: 

• CALPHAD fundamental databases — PrecipiCalc relies on 

CALPHAD-based databases to capture fundamental 

mechanistic features of multicomponent alloys. These 
databases and associated tuning parameters include: 

o Thermo-Calc® [3] compatible thermodynamic databases 
— representing bulk free energy with state variables such 
as composition and temperature; 
o AE — a phase free energy shift in Thermo-Calc to locally 
tune equilibrium phase fractions; 
o DICTRA [4] compatible mobility databases — 

representing atomic mobility and allowing calculation of 
diffusivity when combined with the thermodynamic 
database; 

o D sca i e — diffusivity correction factor used by PrecipiCalc 
to easily rescale the diffusivity matrix for local fitting; 
o Molar Volume — preliminary multicomponent molar 
volume models were developed for both y and y’ under the 
AIM program 

• Material kinetic model parameters: 

o G co h — coherent surface energy, which is the key 
parameter affecting the nucleation barrier when particles 
are small and coherent with the surrounding matrix; 
o G e i — elastic coherency (misfit) energy adds additional 
energy penalty to the precipitation when the particles are 
coherent with the matrix; 

o Gi nC oh — incoherent surface energy captures the increased 
surface energy when the particles lose coherency with the 
matrix; 

o M 0 — prefactor for the interfacial mobility term to 
describe incoherent interfaces, while the corresponding 
activation energy is scaled to that for solvent self diffusion; 
o Rcoh^incoh — the characteristic particle size for 
coherency transition. 

The current PrecipiCalc calibration and validation protocol is 
summarized in Table 2. High-temperature long aging is applied to 
Ni-based superalloys to produce near equilibrium micro structures, 
and matrix/particle compositions are measured by APT and EDS, 
(see next section). The results help down-select a candidate 



thermodynamic database and determine potential simple energy 
corrections AE. Next, diffusion couples between Ni-based 
superalloys and pure Ni with high temperature aging are used to 
determine if there is a need to correct mobility of individual 
components or scale with a simple scaling factor (D sca i e ) to the 
calculated diffusivity from databases. Ultimately, XRD can be 
used to evaluate the importance of misfit and modify the molar 
volume models if needed. 


Table 2 Summary of calibration and validation protocol. 


Experiments 

CALPHAD 

Fundamental 

Databases 

Material Kinetic 
Model Parameters 

Equilibrium Age 
+ APT and EDS 
Compositions 

Thermodynamics, 

AE 


Diffusion Couple 
+ Microanalysis 

Mobility, D sca i e 

SSDTA + APT 


^coh ^ G e i (est.) 

Coarsening Age + 
SEM/TEM for y’ 
size and fraction 


^incoh> M 0 

XRD, TEM for 
misfit 

Molar Volume 

Gel, Rcoh— >incoh 


For material specific model parameters of precipitation kinetics, 
nonisothermal nucleation experiments using SSDTA (see later 
discussion) were identified to determine coherent y’ surface 
energy (a co h) and estimated elastic coherency energy (G e i). Next, 
later stage coarsening experiments are used to determine the 
incoherent y’ surface energy (Gi nco h) and possible interfacial 
mobility term (M 0 ). Finally, TEM study can be used to determine 
the coherency transition size (Rcoh-^incohX and to assist the 
calculation of elastic coherency energy (G e i). 

Thermodynamics 

To assess equilibrium phase fractions and compositions, water 
quenched supersolvus treated samples of ME3, LSHR, Alloy 10 
and RR1000 disc alloys are given lOOOhr treatments at 
temperatures of 1093°C (2000°F), 927°C (1700°F), and 760°C 
(1400°F). Commercial thermodynamic databases compared in this 
study, using Thermo-Calc [3] software, include Thermotech Ni- 
DATA versions 4 to 7 [5], Computherm PanNickel [6], Ni-NIST 
[7] (not including Nb), and TCNI1 [8] (not including Mo, Nb and 
Ta). 

Alloy 10/ME3 Equilibrium Studies 

Comparison in equilibrium phase compositions was undertaken 
with evaluation of the Alloy 10 and ME3 samples quenched from 
the 1093°C lOOOhr age treatment. The high-temperature matrix 
composition of ME3 was determined by both Energy Dispersive 
Spectroscopy (EDS) analysis conducted at NASA GRC and 
Atom-Probe Tomography (APT) using an Imago Local Electrode 
Atom Probe (LEAP) with new larger-FOY LEAP detector at the 
Northwestern University Center for Atom Probe Tomography 
(NUCAPT). In addition, y’ composition was also successfully 
measured by EDS for both alloys. 


A LEAP dataset with 2.2 million atoms and dimensions of 
36x37x60 nm 3 was analyzed for the ME3 sample. No large 
isothermally aged y’ was detected in LEAP (the fine y’ formed 
during quench will be discussed later), and thus the entire LEAP 
dataset represents the matrix during the high temperature age. 
Table 3 summarizes both the quantitative APT and semi- 
quantitative EDS results, and a comparison with the predicted 
high temperature matrix compositions using different 
thermodynamics databases. Based on the root mean square (RMS) 
error of solute content relative to the APT measurement, the Ni- 
DATA 7 database gives the best agreement. In the last column of 
Table 3, an experimental phase fraction of 24.8% was estimated 
by mass balance using the APT measured matrix composition and 
the Ni-DATA 7 predicted y’ composition, and again the Ni- 
DATA 7 prediction provides the best agreement. The comparison 
of predicted versus EDS measured y’ compositions also confirm 
the best agreement from the Ni-DATA 7 database. 

Table 3 Comparison of the predicted y and y’ compositions and 
y’ phase fraction in ME3 at 1093°C (2000°F) with the 
experimental APT y and EDS y and y’ compositions. Ni- 
DATA version 5 (or 7) is abbreviated as Ni 5 (or 7) in 
the table. 




Composition at 1093°C, at% 

% Y 



Ni 

Al 

Cr Co 

Ti 

Mo 

W 

Nb 

Ta 

C 

RMS 

fraction 


APT 

44.9 

6.7 

17.5 21.7 

3.5 

3.4 

1.0 

0.5 

0.5 

0.04 


24.8 


2o 

0.04 

0.02 

0.03 0.03 

0.01 

0.01 

0.01 

0.01 

0.01 

0.002 


1.5 


EDS 

46.1 

6.6 

16.6 22.7 

3.3 

3.0 

1.3* 

0.6 

- 

- 



y 

Ni 5,6 

45.7 

5.6 

19.7 22.1 

2.4 

3.1 

0.7 

0.3 

0.3 

0.01 

0.93 

30.7 


Ni 7 

46.4 

6.0 

18.5 21.8 

2.7 

3.0 

0.7 

0.4 

0.5 

0.01 

0.54 

25.5 


Pan- 

Nickel 

44.3 

5.3 

20.3 23.0 

2.4 

2.9 

0.8 

0.4 

0.5 

0.05 

1.19 

32.3 


NIST-Ni 

44.7 

5.1 

19.9 23.8 

2.4 

3.2 

0.8 


0.2 


1.44 

29.4 


EDS 

58.4 

12.9 

4.1 12.8 

8.9 

0.5 

0.3 

0.6 

1.4 

- 




Ni 5,6 

61.5 

5.7 

2.6 14.4 

7.1 

0.5 

0.9 

2.4 

4.9 

- 

3.09 


Y 

Ni 7 

60.2 

12.0 

2.8 13.0 

8.8 

0.3 

0.3 

1.3 

1.2 

. 

0.63 


Pan- 

Nickel 

61.8 

12.3 

2.2 12.5 

8.2 

0.4 

0.1 

1.2 

1.1 

- 

0.81 



NIST-Ni 

68.3 

6.3 

1.3 9.7 

7.6 

0.3 

0.3 

- 

6.2 

- 

3.27 



* value includes both W and Ta due to peak overlaps in EDS. 


High temperature yand y’ compositions of Alloy 10 (1093°C 
lOOOhr) were successfully measured by EDS, and compared with 
several thermodynamic databases in Table 4. Combined W+Ta 
composition is reported for EDS due to peak overlap. Overall, Ni- 
DATA 7 agrees the best with the semi-quantitative EDS 
measurement. 

A systematic underestimation of Ti in the calculated y matrix was 
also noted in QuesTek’s earlier simulations in the DARPA-AIM 
initiative on commercial Ni-base superalloys [1] and may have to 
be addressed later in this program. 
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Table 4 Comparison of the predicted y and y’ compositions at 
1093°C and EDS measurement for Alloy 10 aged at 
1093°C for 1000 hr. 



Compositions, at% | 


Ni 

Al 

Cr 

Co Ti Mo 

W 

Ta 

Nb 

RMS 


EDS 

53.8 

7.0 

13.4 

17.5 3.2 2.3 

1.6 

1.13 



Ni 5,6 

50.8 

5.74 

17.60 

17.76 2.61 2.30 

2.54 

0.14 

0.50 

1.75 

y 

Ni 7 

51.1 

5.97 

16.97 

17.72 2.76 2.21 

2.48 

0.20 

0.63 

1.49 


MI- 
NI 1ST 

51.3 

5.44 

17.00 

18.60 2.58 2.31 

2.66 

0.10 

- 

1.68 


EDS 

62.3 

13.0 

4.0 

10.9 7.2 0.6 

0.9 

1.3 



Ni 5,6 

62.2 

12.02 

2.82 

11.01 8.42 0.24 

0.93 

0.49 

1.89 

0.80 

Y 

Ni 7 

64.1 

12.11 

2.80 

10.53 8.64 0.23 

0.91 

0.40 

1.73 

0.83 


Ni- 

NIST 

66.8 

13.5 

1.35 

7.54 9.62 0.13 

0.37 

0.68 

- 

1.93 


Ni-Al-Crf-X) (X=Re,W) Model Alloys 


Table 5 Phase composition comparison of APT [9] and the 
database predictions for a Ni-5.2Al-14.2Cr (at%) alloy 
at 600°C. 



600 C 


at% 

APB 

600 C 

% f 



Ni 

Al Cr 

RMS 

(J/m 2 ) 

misfit % 

fraction 



APT 

81.26 

3.13 15.61 






Equilibrium 

matrix 

2a 

0.18 

0.08 0.18 

0.14 





Ni 5,6 

81.40 

3.64 14.96 

0.58 




y 

Ni 7 

TCNI1 

81.34 

3.62 15.04 

0.53 




composition 


81.17 

3.93 14.90 

0.76 






Ni-NIST 

80.60 

5.20 14.20 

1.77 






APT 

76.53 

16.69 6.77 


0.19 

-0.72 

15.4 



2a 

0.50 

0.44 0.30 

0.38 

0.01 


0.4 


Y equilibrium 

Ni 5,6 

74.98 

16.23 8.80 

1.47 

0.19 

-0.6 

12.4 


composition 

Ni 7 

74.97 

16.21 8.82 

1.49 

0.19 

-0.6 

12.4 



TCNI1 

75.59 

14.55 9.57 

2.49 

0.18 

-0.8 

12.3 

Y 


Ni-NIST 

75.96 

12.2511.79 

4.74 

0.15 

-0.9 

0.0 


APT 

72.40 

18.30 9.30 


0.16 

-0.2 



Critical y’- 
nucleus 
composition 

2a 

2.20 

1.80 1.40 

1.61 

0.03 




Ni 5,6 

74.91 

17.37 7.71 

1.30 

0.20 

-0.5 



Ni 7 

74.91 

17.40 7.69 

1.30 

0.20 

-0.5 



TCNI1 

75.36 

16.43 8.21 

1.53 

0.20 

-0.6 




Ni-NIST 

No precipitation 






parameter misfit based on QuesTek's molar volume model. These 
comparisons again support the Ni-DATA 7 database as most 
accurate, and suggest the corresponding error in APB enthalpy 
and misfit is acceptable within current structure/property model 
uncertainty. 

While predicted phase fractions show good agreement with 
experiment for the complex superalloy of Table 3, we note a 
significant discrepancy for the simpler model alloys of Tables 5-8. 
For these alloys a AE shift in the relative free energy of y’ could 
be appropriate for further modeling of precipitation dynamics. 

Table 6 Phase composition comparison of ICP [10] and the 
database predictions for a Ni-10Al-8.5Cr (at%) alloy at 
800°C. 


800 C 

at% 

% Y 

fraction 

Ni Al Cr 

RMS 


ICP 

82.71 8.43 8.86 



y equilibrium 

2a 

0.24 0.16 0.18 

0.17 


matrix 

Ni 5,6 

82.93 8.32 8.75 

0.11 


composition 

Ni 7 

82.93 8.32 8.75 

0.11 



Ni-NIST 

81.90 9.76 8.34 

1.01 



ICP 

76.60 17.41 5.99 


14.6 

y’ equilibrium 

2a 

0.36 0.34 0.12 

0.25 

0.3 

precipitate 

Ni 5,6 

75.52 18.64 5.84 

0.88 

13.9 

composition 

Ni 7 

75.52 18.64 5.84 

0.88 

13.9 


Ni-NIST 

76.47 17.40 6.13 

0.10 

0.0 


Table 7 Phase composition comparison of APT [11] and the 
database predictions for a Ni-Al-Cr-Re alloy at 800°C. 


800 C 

at% 

APB 

(J/m 2 ) 

800 C 
misfit % 

% Y 

fraction 

Ni Cr Al Re 

RMS 



APT 

81.07 10.04 6.74 2.15 






Equilibrium 

motri v 

2a 

0.16 0.04 0.08 0.02 

0.05 




y 

1 1 Id LI IX 

composition 

Ni 7 

80.55 8.94 8.35 2.16 

1.13 






Ni-NIST 

80.22 8.79 8.87 2.12 

1.42 






APT 

76.17 4.97 18.05 0.81 


0.19 

-0.71 

24.8 


Equilibrium 

2a 

0.26 0.10 0.18 0.04 

0.12 

0.00 


0.2 


composition 

Ni7 

75.29 5.68 18.80 0.23 

0.68 

0.20 

-0.68 

15.9 

Y 


Ni-NIST 

75.82 6.13 17.65 0.39 

0.75 


-0.77 

7.2 


composition 

APT 

76.33 5.46 16.92 1.29 


0.17 


20.3 


at 80% 

2a 

0.30 0.16 0.26 0.08 

0.18 

0.00 


0.2 


completion 

Ni7 

75.22 5.56 18.99 0.23 

1.35 

0.20 


12.8 


Further comparison of thermodynamic databases is enabled by 
previous APT and ICP (Inductively Coupled Plasma atomic- 
emission spectroscopy measurement) studies at Northwestern 
University of isothermal y’ precipitation in simple Ni-Al-Cr-X 
model alloys [9-13]. The measured time evolution of particle size 
suggests the precipitation composition trajectory can be modeled 
as an unstable equilibrium with capillarity. Comparison of phase 
compositions with the database predictions for these alloys is 
made in Tables 5-8. For the base Ni-Al-Cr alloy at 600°C, 
measurement of the composition trajectory includes the initial 
critical nucleus composition. As metrics of relative impact of 
phase composition error, Tables 5-8 include not only the RMS 
solute concentration error and y’ phase fraction, but the 
corresponding predicted error in thermodynamically computed y’ 
APB enthalpy [13], and the at-temperature interphase lattice- 
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Table 8 Phase composition comparison of APT [12] and the 
database predictions for a Ni-Al-Cr-W alloy at 800°C. 


800 <C 

at% 

APB 

(J/m 2 ) 

800 <C 
misfit % 

% Y 

fraction 

Ni Cr Al W 

RMS 

y equilibrium 
matrix 
composition 

APT 

2a 

81.31 5.83 11.52 1.34 
0.07 0.04 0.05 0.02 

0.04 




Ni 7 

Ni-NIST 

81.44 6.22 10.53 1.80 

80.48 8.57 8.84 2.11 

0.67 

2.26 




y’ equilibrium 
precipitate 
composition 

APT 

2a 

76.30 17.00 3.91 2.80 
0.08 0.07 0.04 0.03 

0.05 

0.21 

0.00 

-0.49 

37.9 

0.0 

Ni7 

Ni-NIST 

75.39 18.00 4.20 2.42 

75.87 16.61 6.42 1.09 

0.64 

1.77 

0.22 

-0.36 

-0.61 

32.1 

18.7 


Misfit Estimates 


Comparison of relative interphase misfit also provides a useful 
metric of database accuracy for the complex disc alloys. Figure 1 
summarizes the predicted equilibrium at-temperature y/y’ misfit 
for disc alloys, taking into account the presence of equilibrium 
carbides, using phase compositions predicted by the Ni-DATA 7 
database. The alloys fall into 2 groups of relatively high and low 
misfit. Of the 3 rd generation disc alloys, only the RR1000 alloy is 
predicted to maintain the low misfit of the IN100 and Rene88DT 
alloys. In support of this prediction, Figure 2 compares the y’ 
microstructures of supersolvus treated (1193°C/lhr) Alloy 10 and 
RR1000 after long-term aging at 1093°C for 1000 hr. While the 
fine intragranular precipitates in Alloy 10 show the cuboidal 
morphology and ordered arrays consistent with the predicted 
higher misfit, the y’ in RR1000 shows a spheroidal morphology 
and disordered spatial distribution consistent with its predicted 
low misfit. The combination of thermodynamic database and 
misfit model thus appears to give a reasonable ranking of relative 
misfit between alloys. However, the higher misfit values appear 
susprisingly large and XRD studies are planned to calibrate the 
molar volume model. Elastic energy based misfit estimates from 
relative nucleation undercoolings are discussed later. 



Figure 1 Summary of predicted equilibrium at-temperature y/ y’ 
misfit for disc alloys, taking into account the presence 
of equilibrium carbides, using phase compositions 
predicted by the Ni-DATA 7 database and the 
preliminary molar volume model. 


Mobility 

To assess accuracy of available atomic mobility databases for 
combination with thermodynamic databases in order to predict 
multicomponent diffusivities, a linear diffusion “multiple” was 
prepared at NASA GRC. A 5mm (0.2-inch) thick disc of pure Ni 
was diffusion bonded by hot pressing to a 3.2mm (0.125-inch) 
thick disc of ME3 on one side and to a 3.2mm (0.125-inch) thick 
disc of Alloy 10 on the other, for 4 hours at 877°C (1610 °F) in 
vacuum. After the initial bond was formed, the diffusion couple 
was subjected to additional annealing in a horizontal tube furnace 
in an argon atmosphere at two temperatures, 1093°C or 927°C for 
100 or 300 hours. 


DICTRA [4] simulations with Ni-DATA 7 and the NIST mobility 
database [14] of the diffusion multiple at the hot-pressing 
temperature 877°C predicted no significant diffusion. Therefore, 
subsequent simulations ignored the hot-pressing step. 



Figure 2 Comparison of the y’ micro structures of (a) Alloy 10 
and (b) RR1000 (courtesy of L. Birrell at University of 
Cambridge) 


Next, the NIST mobility database employed in the AIM project 
was compared with the Thermo-Calc MOB 1 database in 
combination with various thermodynamic databases. In order to 
compare microprobe analyses of the diffusion multiple with 
DICTRA simulations, the average Matano interface (defined as 
the interface across which equal number of atoms have crossed in 
both directions) was equated to the origin of the calculated 
profiles, and simulations assumed that y’ precipitates act only as 
sink or source of solute for diffusion, i.e. no diffusion through y’ . 
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-300 -200 -100 0 100 200 300 
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Figure 3 Comparison of microanalysis (discrete points) and 
DICTRA (curves) predicted composition trace of (a) 
ME3/Ni diffusion couple aged at 1093°C for 300 hr, 
and (b) Alloy 10/Ni diffusion couples aged at 1093°C 
for 100 hr. 

The agreement with the experimental microanalysis results is best 
using Ni-DATA 7 and NIST mobility databases. However, there 
are two discrepancies: (1) with the exception of Al, the 
experimental data does not confirm the predicted strong 
nonmonotonic profiles near the interface, and (2) the actual 
diffusion of Nb, Ta, Al, and Ti is less than predicted. Using the 
NIST mobility database with minor adjustments to the diffusivity 
prefactors of Cr, Nb and Ti (by factors of 1.4, 0.53 and 0.70, 
respectively), the predicted DICTRA composition profiles agree 
reasonably well with the microanalysis, as shown in two of the 
selected examples in Figure 3. However, comparative simulations 


reveals that these adjustments do not affect the y’ precipitation 
simulation with PrecipiCalc. Hence, we conclude there is no 
correction needed to the NIST mobility database and thus D sca i e 
is set to 1 . 

Furthermore, as seen in the micrograph of Figure 4(b), the width 
of significant grain growth in ME3 due to MC/y’ dissolution also 
agrees well with the DICTRA predicted phase fraction 
distribution in Figure 4(a). 


Phase 

Fraction 


1 52 |Llm 


FCC 


0.1 _ 


(a) 


0.01 


:jvr_c 


0.001 


-400 -300 -200 -100 


100 



Distance (microns) 


Figure 4 Comparison of (a) predicted equilibrium phase 
fractions, and (b) measured phase boundaries with 
adjusted Matano interface, for ME3/Ni 1093°C/100hr 
diffusion couple sample. 


Nonisothermal Precipitation 

Nucleation Onset Temperatures and SSDTA 

As an efficient experimental check of nucleation behavior, 
transformation onset temperatures were measured in rapidly 
cooled arc melted samples of ME3 and Alloy 10 using a Single 
Sensor Differential Thermal Analysis (SSDTA) technique 
developed at Ohio State University [15]. From nucleation theory, 
the most important material parameters affecting the nucleation 
rate are driving force (which is determined by the bulk 
thermodynamics), surface energy and elastic coherency energy. In 
this section, utilizing the Ni-DATA 7 and NIST mobility 
database, we present the results of calibrating the surface energy, 
while initially ignoring the coherency elastic energy. 

The transformation onset temperatures (T onset ) from the SSDTA 
experiments were determined from SSDTA data processing 
software, which takes the measured temperature versus time, fits 
the baseline prior to the transformation with a choice of function 
suitable for describing cooling, and then determines the departure 
point of measured data from the fitted function. The calibration of 
the SSDTA thermocouples has shown an accuracy of ±4.5°C [16] 
in the relevant temperature range. To determine the onset 
temperature from PrecipiCalc simulation with physically 


927 



equivalent interpretation as closely as possible to the SSDTA 
results, the following summarizes our procedure. 

• Process SSDTA measured temperature profiles to remove the 
latent heat contribution; 

• Calculate material compositions relevant to the y’ 
precipitation from matrix y phase — the high temperature 
phases (borides, carbides and undissolved y’) are removed 
with equilibrium calculations at highest SSDTA measured 
temperatures using the Ni-DATA 7 database. 

• Perform PrecipiCalc simulations using Ni-DATA 7 and NIST 
mobility databases, with estimated surface energy; 

• Collect time (or temperature) evolution results of y’ volume 
fraction, and compositions of matrix and y’, compute the time 
(or temperature) evolution of molar enthalpy; 

• Compute temperature derivatives of the molar enthalpy with 
respect to temperature, dH/dT; 

• Determine the transformation onset temperature where dH/dT 
changes by more than 10% of average dH/dT values, with 
decreasing temperature. 

We note the difference of calculated T onset values with the above 
enthalpy change approach and a simple approach of using a small 
y’ volume fraction of 0.1% is within the SSDTA accuracy. 



Average Cooling Rate (C/sec) 

Between 1160 and 1125C 

Figure 5 Summary of SSDTA results and PrecipiCalc (PpC) 
predictions with calibrated coherent interfacial energy 
o CO h (values in J/m 2 ). The data are plotted with respect 
to average cooling rate from measured temperature 
profiles. 


Figure 5 summarizes the interfacial energy calibrated PrecipiCalc 
results and comparison to the SSDTA measured T onset , for all 4 
alloy samples. We varied the surface energy within 0.0225 to 
0.0325 J/m 2 , and Figure 5 shows the best fit surface energy and 
T 0 nset results with minimum RMS of the difference between 
prediction and measurement. Except LSHR (which will be further 
discussed later), the optimized interfacial energies of the 
remaining three alloys are similar and agrees with the value of 
interfacial energy, 0.022-0.023 J/m 2 , reported by Sudbrack et.al. 
[9] for the Ni-Al-Cr ternary model alloys. 

The experimental uncertainty of T onset associated with SSDTA is 
represented by error bars in Figure 5. The PrecipiCalc predictions 
with optimized interfacial energy agree with SSDTA results well, 
and mostly fall within the experimental uncertainty. Note that the 
Tonset temperatures, both measured and predicted, do not 
monotonically decrease with the average cooling rate between 
1 125 to 1160°C. If the temperature followed linear cooling, T onse t 
should have decreased monotonically with higher cooling rate. 
The actual local temperature profile clearly affected the measured 
Tonset’ an d PrecipiCalc predictions capture the non-monotonic 
behavior very well. 

Due to a furnace temperature limitation of the SSDTA 
experiment, Alloy 10, having the highest predicted y’ solvus 
temperature, contains undissolved large y’ as confirmed by SEM 
(Figure 6). The observed amount is in good agreement with 
predicted equilibrium y’ fraction at the highest sample 
temperature in the furnace. The calculated T onS et °f Alloy 10 
presented in Figure 5 was based on the assumption that 
undissolved y’ does not grow during the quench, allowing us to 
remove the large y’ completely in the PrecipiCalc simulations. To 
test this assumption, we performed an Alloy 10 simulation 
allowing primary y’ growth during the quench, which resulted in a 
small reduction of nucleation onset temperature by only 1.6°C. 
Hence, we conclude the assumption to remove large undissolved 
y’ in SSDTA Alloy 10 quench simulation is reasonable. 



Figure 6 NASA SEM observation of an AlloylO SSDTA sample. 

y’ volume fraction was estimated to be around 9.5%, 
which is compatible with Ni-DATA 7 equilibrium 
prediction of 8.4% y’ phase. 

The relative ranking of best-fitting interfacial energies in Figure 5 
ranks with the relative experimental nucleation undercoolings. 
This in turn ranks with the relative predicted misfit (RR1000 
lowest, LSHR highest). Incorporation of a coherency elastic 
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energy G e j in further modeling will likely reduce the range of the 
coherent a values. Using the surface energy of 0.0225 J/m 2 , found 
for RR1000 and assumed to correspond to zero misfit, we reran 
PrecipiCalc for ME3, LSHR and Alloy 10 for SSDTA 
simulations with the same surface energy and then estimated the 
coherency elastic energy G e j by fitting to the measured nucleation 
onset temperature. Table 9 summarizes the results, and also the 
corresponding at-temperature misfits based on the theoretical 
elastic strain energy of inclusions [17]. According to this analysis, 
LSHR corresponds to the highest misfit, followed by ME3 and 
Alloy 10, in a similar ranking but at about 1/3 of the misfit values 
of the model prediction in Figure 1 . These misfit estimates will be 
validated against ongoing XRD experiments. In subsequent 
PrecipiCalc simulation to be discussed here, the surface energies 
in Figure 5 are used. 

Table 9 Elastic energy based misfit estimate from SSDTA 
nucleation onset temperatures and RR1000 surface 
energy. 



Estimated G e i 

Misfit 


(J/mole) 

(%) 

RR1000 

0 

0 

Alloy 10 

30 

0.5 

ME3 

50 

0.61 

LSHR 

80 

0.8 


Y Micro structure in SSDTA Samples 


The well-defined thermal history of the SSDTA samples also 
provides an excellent opportunity for experimental validation of 
PrecipiCalc simulations of nonisothermal precipitation. Using 
LEAP APT, we analyzed the ME3 Ar gas-quenched sample 
(ME3-lAr) from the SSDTA experiment. In total, 32 million 
atoms were analyzed. The overall bulk composition is very close 
to the expected composition. Larger y precipitates, about 40- 
60nm in diameter, were observed together with ultrafine y’ 
precipitates, about 5nm or less in diameter. 

High-resolution SEM (HR SEM) analysis was also conducted on 
the same ME3-lAr SSDTA sample, as well as an LSHR Ar gas- 
quenched sample (LSHR-2Ar). Consistent results between the 
APT and SEM for ME3-lAr are shown in Figure 7. As shown by 
the higher magnification views at the bottom of Figure 7, fine y’ 
formed at lower temperature during SSDTA cooling was also 
observed in both ME3-lAr LEAP APT and HR SEM. 

Table 10 shows a summary of measured larger y’ size and also 
PrecipiCalc predictions using the measured SSDTA temperature 
profiles. PrecipiCalc provides reasonable estimates with the 
approximations employed thusfar. Figure 8 shows reasonable 
agreement in particle size histograms between HR SEM and 
simulation for the LSHR-2Ar sample. In RR1000-3He (He gas- 
quenched) LEAP study, fraction of 35% was determined for 
largest y\ which agrees well with PrecipiCalc’ s 37.6%. 

As in the case of the model alloy of Table 8, the APT 
measurements in the ME3-lAr sample show W partitioning in the 
direction of thermodynamic prediction for the larger y’ 
population. However, the W profiles observed for the ultrafine 
population suggest solute trapping of the slowest diffusing 


components may be occurring during precipitation at the lowest 
temperatures. 



Figure 7 Characterization of y’ micro structure in ME3-lAr 
SSDTA sample. Left (a,c) figures from LEAP APT and 
right (b,d) figures from HR-SEM are compared at two 
different length scale — (a) and (b) at smaller 
magnification, showing larger y’, and (c) and (d) at 
slightly higher magnification showing very fine y’ . 



y Diameters (nm) 


Figure 8 Comparison of particle size histogram of PrecipiCalc 
(3D results) and HR SEM with an observation of 30 y’ 
particles in an LSHR-2Ar SSDTA sample. 
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Table 10 Comparison of the largest y’ mean diameters between 
two SSDTA samples and PrecipiCalc simulation using 
the coherent interfacial energies discussed previously. 



LEAP 

APT 

HR SEM 

PrecipiCalc 

Predictions 

ME3-1 Ar 

40-60 nm 

53.8 nm 

41 nm 

LSHR-2Ar 

- 

50.1 nm 

45.6 nm 

RR1000-3He 

20-45 nm 

31.8 nm 

20.5 nm 


y’ Micro structure in ME3 Water Quenched Samples 

LEAP APT microanalysis was conducted on ME3 samples water 
quenched from both subsolvus treatment at 1093°C (lOOOhr) and 
supersolvus treatment at 1200°C (lhr). For the supersolvus treated 
sample, a LEAP dataset containing 8.7 million atoms within a box 
of 128x55x54 nm 3 was collected. The LEAP data analysis shows 
nanoscale y’ precipitates of about 15~30nm in diameter, which is 
in good agreement with our PrecipiCalc prediction of 30nm with 
an estimated cooling temperature profile. TEM work at NASA 
GRC indicates a larger mean diameter of y’ at 44nm. 

For the 1093°C subsolvus treated ME3, the LEAP experiment 
used to obtain 1093°C matrix composition (Table 3) was further 
analyzed for fine quench y\ An estimated 10-20nm effective 
diameter from LEAP APT is in reasonable agreement with the 
predicted 29nm diameter from PrecipiCalc simulation with an 
estimated cooling temperature profile. A proximity histogram of 
the subsolvus LEAP sample summarizes the measured 
partitioning of solutes between the y and y’ phases in Figure 9. As 
summarized in Table 11 for both subsolvus and supersolvus 
treated samples, the partitioning is again in the direction of 
thermodynamic prediction, with the notable exception of the flat 
W profile offering further evidence that the slowest diffusing 
component exhibits solute trapping at extreme cooling conditions. 
Further study will address whether solute trapping need be 
considered for industrially-relevant processing conditions. 



Figure 9 y and fine quench y’ compositions (in atomic %) 
analysis of LEAP experiment on ME3 1 093 °C/ lOOOhr 
isothermal hold sample. 


Table 1 1 Comparison of measured nanoscale y’ particle and 
matrix compositions in ME3 with corresponding 
PrecipiCalc simulations. 


ME3 

Compositions, at% 

Ni Al Cr Co Ti Mo W Nb Ta 

RMS 

Subsolvus 

42.5 4.9 21.3 23.7 2.3 3.6 1.0 0.4 0.3 

35.0 2.82 27.3 27.7 0.51 5.28 1.30 0.07 0.08 

2.80 

Y 

APT 

PpC 

Y 

ATP 

PpC 

55.6 9.3 7.3 17.5 6.1 2.1 0.8 0.4 0.8 

61.3 12.8 2.1 12.2 8.2 0.42 0.52 1.2 1.2 

3.07 

Supersolvus 

40.5 4.8 22.5 23.8 2.0 4.2 0.8 0.6 0.5 

39.2 3.00 24.6 27.5 0.74 3.94 0.87 0.08 0.13 

1.71 

y 

APT 

PpC 

Y 

APT 

PpC 

56.5 12.0 3.5 13.6 8.9 1.9 0.5 1.0 1.6 

62.2 12.4 1.94 11.9 8.50 0.32 0.34 1.0 1.3 

1.01 


y’ Micro structure in Furnace Cooled Samples 

Supersolvus treated ME3, LSHR and Alloy 10 samples followed 
by instrumented furnace cooling were characterized under SEM. 
The cooling rate between 1160 and 1125 °C is around 0.54C/sec, 
which is about two orders of magnitude slower than the SSDTA 
and water quench experiments. These experiments offer growth- 
dominated conditions for calibration of interfacial mobility. 
PrecipiCalc simulations were conducted using the coherent 
surface energies determined from the SSDTA study and the 
individually measured temperature cooling profiles for each 
samples. The PrecipiCalc 3D results converted to 2D are in 
reasonable agreement with the SEM results, as shown in Table 12. 
Unlike the finer y’ particles formed during SSDTA fast cooling 
(see Figure 8) whcih were etched to reveal 3D sizes in SEM, y’ 
particles here are much larger and should be compared to the 2D 
results. As shown in Figure 10, we observe reasonable agreement 
between simulation and experimental size distribution for these 
growing supercritical y’ particles, showing sharper size 
distribution than that from coarsening theory. Note the truncation 
of the calculated 2D distribution is lOnm, which is about the 
resolution of SEM and image analysis. 

These PrecipiCalc simulations were carried out with coherent y/y’ 
surface energy and without interfacial dissipation (M 0 term in 
Table 2). Though these y’ particles are significantly larger than 
those found in SSDTA and water quench experiments, there is no 
indication that they are incoherent with the matrix. However, it is 
possible that they are semi-incoherent and are in the process of 
losing coherency. 

Table 12 Comparison of the largest y’ mean sizes (diameters) 
between furnace cooled samples and PrecipiCalc 
simulations. 



HR SEM 

PrecipiCalc Results 
converted to 2D 

ME3 

430 nm 

348 nm 

LSHR 

408 nm 

402 nm 

Alloy 10 

376 nm 

396 nm 


930 




4 


Q 3 

CO 

CL 


CO 

E 

o 

2 1 

0 

0 0.5 1 1.5 2 

R/<R> 

Figure 10 Comparison shows good agreement between 
experimental SEM (bars) and PrecipiCalc simulated 
(converted to 2D, solid line) and normalized particle 
size distributions for LSHR furnace cooled sample. 2D 
coarsening size distribution is shown in a dash line. 



growth leads to the sharper size distribution of Figure 10. The 
good agreement with average size in Table 12 suggests that 
interfacial dissipation can be negliected for particle size up to the 
500nm scale. 


Coarsening Experiments 

ME3, LSHR and Alloy 10 samples were supersolvus treated at 
1193°C for 1 hr, water quenched, and then aged at different 
conditions, followed by a water quench. The aging conditions 
selected are 1093°C for 20 and 1000 hr, 927°C for 1000 hr, and 
760°C for 1000 hr. The samples were then analyzed under TEM 
or SEM, depending on the size scale of y’ particles. For higher 
temperature aging, both intragranular y’ and grain boundary y’ 
particles are monitored as represented in in Figure 12. Similar to 
Alloy 10 in Figure 2a, the intragranular y’ in ME3 shows an 
ordered alignment, indicating coherency based elastic interaction. 
A consistent image processing technique was used to identify the 
grain boundary y’ and to obtain the sizes and fractions of both 
types of y’ particles for all ME3, LSHR and Alloy 10 samples. All 
of these heat treatment and experimental analysis were conducted 
at NASA GRC. 




Figure 1 1 Comparison of PrecipiCalc simulated y’ (a) mean 
radius (<R>), (b) growth rate and (c) mean radius over 
critical radius (R c ) ratio with respect to temperatures 
for both LSHR furnace cooled (FC) and SSDTA (2Ar). 
In this temperature range, only the first y’ distribution 
is formed in these two calculations. 

Figure 1 1 compares the PrecipiCalc simulated y’ size and growth 
rate evolution of fast cooled SSDTA and slow furnace cooled 
LSHR samples. In these two very distinct cooling rates, y’ 
nucleates at similar sizes within a narrow temperature window 
(1 140-1 142°C). In the slower furnace cooled case (LSHR-FC), 
nucleated y’ particles can grow to much higher multiples of 
critical size (Figure lie), despite the slower growth rate (due to 
lower supersaturation in slower cooling) compared to the fast cool 
at temperatures below 1130°C (Figure lib). This dominance of 



i i i i i i i i i i i 


1764-ME-SW 20.0kV 12.0mm xl.OOk SE(L) 5/14/2007 30.0um 


Figure 12 SEM image (processed to show different features) of 
ME3 supersolvus treated and then aged at 1093°C for 
lOOOhr. Intragranular and grain boundary y’ are 
represented as black and grey, respectively. 

Table 13 summarizes the measured y’ sizes in these heat treatment 
conditions. The grain boundary y’ particles are typically larger 
than the intragranular y’ particles by a factor of 1.5 to 2. Between 
ME3, LSHR and Alloy 10, there are no consistent trends in the 
relative y’ sizes. We note that the TEM and SEM measured sizes 
correspond to the mean 3D and 2D sizes, respectively. 

The PrecipiCalc modeling of coherent intragranular y’ in these 
isothermally aged coarsening experiments is ongoing. Program 
plans include the modeling of incoherent grain boundary y’ using 
PrecipiCalc ' s heterogeneous grain boundary nucleation model. 
The grain boundary y’ growth simulations will include effects of 
grain boundary diffusion and possibly interfacial mobility to 
account for incoherent interfaces. 
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Table 13 Summary of experimental measured f mean equivalent 
diameters (jam) of isothermally aged samples. 
Intragranular (G) f and grain boundary (GB) y\ when 
observed, are reported separately. 


Aging 

Conditions 

Exp 

ME3 

LSHR 

Alloy 10 

G 

GB 

G 

GB 

G 

GB 

Before Aging 

TEM 

0.044 

- 

0.028 

- 

0.044 

- 

1093C-20h 

SEM 

0.510 

0.810 

0.504 

0.984 

0.546 

1.084 

1093C-1000h 

SEM 

1.380 

2.840 

1.420 

2.996 

0.894 

2.618 

927C-1 OOOh 

SEM 

0.584 

0.758 

0.598 

0.954 

0.670 

1.224 

760C-1 OOOh 

TEM 

0.058 

- 

0.062 

- 

0.062 

- 


Conclusion and Summary 

Measured y- y’ phase compositions in equilibrated 3 rd generation 
disc alloys and reported phase compositions in simple Ni-Cr-Al-X 
model alloys show best agreement with the Thermotech Ni-Data 
version 7 thermodynamic database. In the content of 
multicomponent disc alloys, the assessed sensitivity of predicted 
phase fractions, APB enthalpy and interphase misfit to measured 
composition error is judged as acceptable within current 
structure/property model uncertainty. Prediction of a relatively 
high misfit in the high-refractory 3 rd generation alloys is 
supported by micro structural evidence, but precise predictions 
will require further misfit model calibration. Comparison of 
measured composition profiles in diffusion couples of ME3 and 
Alloy 10 with pure Ni against DICTRA multicomponent diffusion 
simulations indicates the combination of the Ni-DATA version 7 
thermodynamic database with the NIST mobility database gives 
sufficient accuracy to support PrecipiCalc simulation of 
diffusional precipitation without diffusivity modification. The 
SSDTA measurement of critical nucleation undercoolings defines 
reasonable values of coherent interfacial energies. Incorporation 
in PrecipiCalc simulation of nonisothermal precipitation for both 
supersolvus and subsolvus treated materials gives good agreement 
with observed particle sizes and compositions. APT microanalysis 
shows some evidence of solute trapping of W, the slowest 
diffusing component, in fine precipitates forming at extreme 
cooling rates. Further modeling of measured long-term coarsening 
behavior requires refinement of coherency elastic energy 
estimates and treatment of the parallel evolution of incoherent 
precipitation at grain boundaries. The accuracy of predictions 
validated so far indicates that available fundamental databases and 
precipitation models offer sufficient fidelity for effective 
application of the AIM methodology to the accelerated 
optimization of hybrid thermal processing of 3 rd generation 
aeroturbine disc alloys. 
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